Breeding and Solitary Wave Behavior of Dunes 
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Beautiful dune patterns can be found in deserts and along coasts due to the instability of a plain 
sheet of sand under the action of the wind. Barchan dunes are highly mobile aeolian dunes found in 
areas of low sand availability and unidirectional wind fields. Up to now modelization mainly focussed 
on single dunes or dune patterns without regarding the mechanisms of dune interactions. We study 
the case when a small dune bumps into a bigger one. Recently Schwammle & Herrmann {Solitary 
wave behavior of sand dunes. Nature 426, 610 (2003)) and Katsuki et al. (Collision dynamics of 
two Barchans dunes simulated by a simple model, cond-mat 04033121 ('2004')') have shown that under 
certain circumstances dunes can behave like solitary waves. This means that they can "cross" each 
other which has been questioned by many researchers before. In other cases we observe coalescence, 
i.e. both dune merge into one, breeding, i.e. the creation of three baby dunes at the center and 
horns of a Barchan, or budding, i.e. the small dune, after "crossing" the big one, is unstable and 
splits into two new dunes. 
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PACS numbers: 45.70.Qj,45.70.Vn,89.20.-a 

I. INTRODUCTION 

We observe many different dune patterns in nature, as 
for example longitudinal, transverse, star and Barchan 
dunes. In regions where the wind blows mainly from the 
same direction sand availability determines the dune pat- 
tern. At high sand disposal transverse dunes dominate 
in the fields. They seem to be translationally invariant 
so that the lateral sand flux can be neglected. 'When less 
sand is available, Barchan dunes appear (Fig.^. These 
are highly mobile having the form of a crescent moon. 
Their velocity can reach up to several tenths of meters per 
year and is proportional to the reciprocal height, mean- 
ing that smaller dunes are faster than large ones. The 
surface of a Barchan can be divided into different sec- 
tions: the windward side, the slip face after the brink and 
the horns from which sand can leave the dune. Barchan 
dunes of different sizes are not perfectly shape invari- 
ant and there exists a minimal height of 2-3 meters be- 
low which they are not stable. These observations result 
from field measurements that have been made over the 
last decades [ll |2, 13, U, S la LD| • Still many questions on 
dune dynamics remain open. Due to the large time scales 
involved in dune formation, full evolution of a dune is 
difficult to assess through measurements. Attempts have 
been made to get more insight through numerical cal- 
culations. Recently, several numerical models have been 



proposed in order to explain dune morphology and for- 
mation |iliE3,ElEiIIiEiEiEiEilmiia,;^EII 

They have to deal with the calculation of the turbulent 
wind field, the saltation sand flux over the windward side 
and the avalanches going down the slip face. Up to now 
modelization mainly focussed on single dunes or dune 
patterns without regarding the mechanisms of dune in- 
teractions. 

Recently, Besler found small Barchans at the down- 
wind side of big ones and concluded that Barchan dunes 
could behave like solitons j22i, 23] . This means that they 
would behave like solutions of non-linear equations, for 
example those describing waves in shallow water, which 
propa gate through each other without changing their 
shape [2J|. As an example see Fig. ^ In (a) and (b) a 
small Barchan is apparently ejected from the main dune, 
whereas, in (c) small dunes emerge from the horns. Note 
the similarity with the snapshots of a collision simula- 
tion depicted in (b). Similar occurrences can be found in 
experiments with sub-aqueous Barchans |27|. Neverthe- 
less, most researchers believe that if a small Barchan hits 
a bigger one, it will be completely absorbed. This is mo- 
tivated by the fact that a sand formation cannot cross the 
slip face of a dune without being destroyed. Therefore the 
description of Barchan dunes as solitons has found very 
little support up to now, until Schwammle et al. found 
that dunes can behave as solitary waves under certain 




FIG. 1: Examples of the breeding process in two Barchan fields, Namibia (a) and Peru (b), (c), and in numerical simulations 

(c). 



conditions 25]. They show that, due to mass exchange, 
a big Barchan colhding with a smaller one placed be- 
hind, may decrease its height until it becomes smaller, 
and therefore, faster, than the previous one, and leaves. 
Meanwhile the initially smaller dune increases its height 
becoming bigger and slower, in such a way that it looks 
as if the smaller dune crosses the big one. This situation 
was referred as solitary wave behavior. Katsuki et al. 
[26j also have obtained solitary wave behavior for coaxial 
and offset collisions of two sub-aqueous Barchans. 

In the following we use a minimal model for dunes to 
develop the morphological phase diagram for the coaxial 
collision between two Barchans. We also characterize the 
transition between the different regimes using two order 
parameters and show a possible size selection mechanism 
in Barchan fields. 



II. MODEL 

Our model [la,lji?i 2l'| consists of three coupled equa- 
tions of motion calculating the shear stress of the wind 
field, the sand flux , the avalanche flux and the resulting 
change of the topography using mass conservation. The 
shear stress of the wind is obtained from the perturbation 
of the air flow over a smooth hill using the well known 
logarithmic velocity profile of the atmospheric boundary 
layer ^l|. The shear velocity describes the strength of 
the wind. At the windward side of a Barchan dune the 
sand grains are transported by the so called saltation 
mechanism [l|. Below the big vortex of the wind field 
behind the brink the air cannot maintain sand transport. 
So the grains are deposited behind the brink in the lee 
zone of the Barchan until the surface reaches the angle of 



repose on which the sand flows down the slip face through 
avalanches. 

The simulations are carried out with a completely uni- 
directional and constant wind source. In every iteration 
the horizontal shear stress r of the wind, the saltation 
flux q and the flux due to avalanches are calculated. The 
time scale of these processes is much shorter than the 
time scale of changes in the dune surface so that they are 
treated to be instantaneous. We perform all simulations 
using open boundary conditions with a constant influx. 
In the following the different steps at every iteration are 
explained. 

a. The air shear stress t on the ground: The shear 
stress is computed according to an analytical work de- 
scribing the perturbation of the ground shear stress by a 
low hill or dune 28] . This perturbation is given by: 
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where a = ^i L k^ zi^jl. 

Here, Kq and K\ are modifled Bessel functions, k^ 
and ky are the components of the wave vector fe, the 
coordinates in Fourier space, t^ and Ty are the Fourier- 
transformed components of the shear stress perturbation 
in wind direction and in transverse direction, h^ is the 
Fourier transform of the height profile, U is the verti- 
cal velocity profile which is suitably nondimensionalised, 
I the depth of the inner layer of the flow, and zq the 
roughness length which takes into account saltation. L 
is a typical length scale of the hill or dune and is given 
by a quarter of the mean wavelength of the Fourier rep- 
resentation of the height profile. In order to take into 



account the shear stress strong reduction in the lee side 
of the hill, where the linear models overestimate it [29j, 
we multiplied the imaginary part of fj. by the fenomeno- 
logical constant 1.5. 

It has to be taken into account that the flow separates 
at the brink of the slip face of a dune. This is done 
by assuming an idealized "separation bubble" , a region 
inside which there is no flow and outside of which the air 
flows as over a shape of the combined dune and bubble. 
Each slice in wind direction of the bubble is modelled by 
a third-order polynomial so that in the case of a barchan 
the region between the horns is inside the bubble. The 
height profile hg in (Q is the profile of the dune including 
the separation bubble. 

b. The saltation flux q: From the shear stress, the 
modification of the air fiow due to the presence of saltat- 
ing grains is accounted for. This results in an effective 
wind velocity driving the grains. 
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where m* = i/r/pair 

The shear stress r results from (^ through t = tq + 
|to| T, where tq is the shear stress over a flat plane, zq 
is the roughness length of the sand excluding the effect 
of saltation, and k w 0.4 is von Karman's constant. Zm, 
the mean saltation height, and zi are parameters of the 
model. 

The next step is the computation of the typical velocity 
of the saltating grains. It is determined by the balance 
between the drag force acting on the grains, the loss of 
momentum when they splash on the ground, and the 
downhill force. It is computed by solving the following 
quadratic vector equation numerically: 
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where Voff — Wcff ^77 ■ Cd is the drag coefficient of a grain, 
engrain sud /Og,ai„ arc the diameter and the density of the 
grains, a is a model parameter. 

The sand flux due to saltation is then obtained by nu- 
merically solving the transport equation: 
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Here g is the gravity acceleration, and a, j3 and 7 are 
model parameters taken from |82| . 

c. The time evolution of the surface: When the sand 
flux has been calculated, the height proflle is updated 
according to the mass conservation: 
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Psand 

d. Avalanches: In the last step, avalanches are sim- 
ulated where necessary. If the slope of the sand surface 
exceeds the static angle of repose, sand is redistributed 
according to the sand flux: 



Qavai = ^ (tanh |V/i| — tanh(tan6'dyn) 
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The surface is repeatedly changed according to ® using 
this flux, until the maximum slope is below the dynamic 
angle of repose, ^dyn- The hyperbolic tangent function 
serve only to improve convergence. 

All these steps are repeated iteratively to simulate the 
evolution of the shape. 

III. RESULTS 

We performed calculations by numerically solving the 
set of equations initially placing a big Barchan (vol- 
ume Vh) downwind of a smaller one (volume Vh). The 
strength of the wind blowing into the system is fixed to 
a shear velocity of 0.5 m/s. The influx is 0.001 kg/ms, 
equal to the big Barchan equilibrium outflux. In order to 
take into account the lack of scale invariance of Barchans, 
we repeat the simulations for two different sizes of the big 
Barchan, with initial volumes Vh- 6 and 70 x lO^m^. The 
same general picture was observed. The smaller Barchan 
at some point bumps into the larger one. This leads to 
a hybrid state where the two dunes melt into a complex 
pattern. Four different situations can be observed: coa- 
lescence (Fig. |2Ji), breeding (Fig. |23 and ^), budding 
(Fig. [21i) and solitary wave behavior (Fig. [2^) depending 
only on the relative sizes of the two dunes. Thus, we 
chose as control parameter the relative volume between 
the two dunes Vh/Vh- 

The evolution of the hybrid state can be understood as 
the result of a competition between two processes. The 
first one is the overlapping of both dunes at the beginning 
of the collision that eventually can lead to coalescence 
(Fig. 13 upper part) . The second one is the effective mass 
exchange between the dunes due to the changes induced 
to the wind shear stress due to the approaching of both 
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FIG. 2: Different situations during the collision of two Barchans for (Vh/Vf/), = 0.06 (a), 0.08 (b), 0.12 (c), 0.17 (d) and 0.3 
(e) using open boundary condition. Coalescence (a), breeding (b) and (c), budding (d) and solitary wave behaviour (e) take 
place. The time (t) is in month. The initial volume and height of the big Barchan is 6 x lO^in^ and 5 m height, whereas the 
heights of the smaller Barchan are 1.8m (a), 1.9m (b), 2.2m (c), 2.6m (d) and 3.1 m (e). 



dunes. In the hybrid state the wind shear stress over the 
windward side of the bigger dune is reduced and thus, 
crest erosion is enhanced. Besides, the wind shear stress 
over the lee side of the dune smaller is also reduced but 
enhancing crest deposition. Thus, the dune smaller may 
gain enough sand to become bigger than the one in front 
and therefore also becomes slower. In this way, the before 
bigger dune can become the smaller one and its velocity 
sufficiently large to leave the hybrid state. In this case 
the dunes separate again (Fig. O bottom) . 

The collision process is crucially affected by the sepa- 
ration bubble, i.e. the region after the brink and between 
the horns at which flow separation occurs (see item a. in 
the Model section) . After the separation at the brink, the 



flow streamlines reattach smoothly near the line segment 
whose end points are the horns. There, sand transport 
continues again. However, inside the separation bubble 
the flow is strongly reduced and, for simplicity, we set the 
flux to zero. Hence, the upwind dune will absorb that 
part of the downwind dune inside its separation bubble 
(Fig. 13;, |3iand|l;). 

For small relative volume {Vh/Vn < 0.07) both dunes 
coalesce to a single one. In this case the relative velocity 
is high and hence the overlapping is faster than their mass 
exchange. Small dunes have a short slip face which dis- 
appears while climbing up the bigger one. This reduces 
the mass exchange and leads to a complete absorption 
of very small Barchans (Fig. 13^). For larger Vii/Vh the 




FIG. 3: Time evolution of the central slice of the 'breeding' 
collision represented in the figure ^p. Note initially the small 
dune climbing on the biger one, and finally their mass ex- 
change that leads to their separation. 



slip face survives for longer time, mass exchange becomes 
relevant, and a small barchan is ejected from the central 
part of the dune (Fig. |5J|. The perturbation of the big 
dune shape, due to the overlapping of the small dune be- 
hind, also propagates over the horns since there is no slip 
face. At the end of each horn a small dune is ejected. 
This phenomenon of triple ejection we call 'breeding'. 
Fig. [2l3 shows the snapshots. Note the qualitative simi- 
larity with the Barchan field shown in Figure ^ 

As the relative volume increases {Vh/Vn > 0.14), a 
smaller relative velocity favours the mass exchange and 
reduces the overlapping process, leading to the complete 
separation of both dunes. Nevertheless, the leaving dune 
lacks the central part of its windward side and cannot 
reach the stable Barchan shape. Therefore, it splits 
into two new dunes, a phenomenon we call 'budding' 
(Fig. 121) . A similar phenomena was reported in experi- 
ments with sub-aqueous Barchans |27l| . 

Between the 'breeding' and 'budding', a transition oc- 
curs at which the 'breeding' three ejected dunes are con- 
nected forming a single dune that, afterwards, splits into 
three again (Fig. O;) . For a higher relative volume the 
central ejected dune decrease its size until disappear at 
'budding'. We consider this transition also as 'breeding'. 

When Vh/Vh is greater than 0.25, the instability of 
the dune in front disappears and both dunes develop to 
Barchans. Then, we observe solitary wave behavior as 
shown in Figure ^. In that case the dunes move with 
similar velocity and the mass exchange is the main pro- 



cess of the evolution of the hybrid state. The overlap- 
ping of both dunes is very small now and the emerging 
dune looses merely a small fraction of its tail. Effectively 
it looks as if the smaller dune just crosses the bigger 
one while in reality due to mass exchange the two heaps 
barely touches each other. 

The morphological phase diagram crucially influences 
the evolution of the number of dunes in a Barchan field. 
In the coalescence region {Vh/Vn < 0.07) the number 
of dunes decreases by one, whereas three dunes in the 
breeding (0.07 < Vh/Vh < 0.14) and two dunes in the 
budding region (0.14 < Vh/Vn < 0.25) appear, and the 
number of dunes increases by two and one respectively. 
Finally, in the solitary waves region (0.24 < Vh/Vn) a 
small dune seems to cross the bigger one and the number 
of them remains constant. 
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FIG. 4: Relation between the relative volume of the dunes 
before ((Vh/Vn)*) and after the collision ((Vh/Vw)/). The 
regions in the morphological phase diagram are showed: co- 
alescence (c), breeding (b), budding (bu) and solitary waves 
(s). For coalescence only one dune emerges and thus, the rel- 
ative volume is zero, whereas for breeding, budding and soli- 
tary waves the relative volume seems to follow a exp(-l/x) 
law (full line). 

In order to analyze the transition type between any 
two regimes, we introduce two order parameters, the 
relative volume after the collision ((V/i/Vtr)/) for the 
'coalescence-breeding' transition, and the inverse of time 
during which the emerging dune does not split, for the 
'budding-solitary waves' transition. Note that the for- 
mer parameter is zero for coalescence, where no dunes 
emerge, and the later one is zero for solitary wave behav- 
ior, where the leaving dune evolves into a Barchan shape 
without splitting. 

Figure 121 shows the relative volume after the collision 



as a function of the initial one, approximately related by 
the formula: 
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where rj, = 0.03. This suggests an essential singularity in 
the 'coalescence-breeding' transition. It implies that the 
volume of the leaving dune in the budding regime does 
not have a minimum value, thus, after a Barchan it is 
possible to find small dunes of any height. 
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FIG. 5: Inverse of the time during which the emerging dune 
does not split, as a function of the initial relative volume. 
By definition, this time is infinitely long for the solitary wave 
behavior. Note the continuity at the transition from budding 
to solitary waves, suggesting a second order transition. 



Thus, we cannot define in a Barchan field a characteristic 
length for the budding process. 

It is interesting to note that the relative volume in- 
creases after a solitary wave collision (Fig. 01 . The col- 
lision process redistributes the initial mass making both 
dunes more similar, giving rise to a size selection mech- 
anism in Barchan fields jl3|. However, due to the per- 
manent exchange of sand between the dunes, their sizes 
are not constant after the collision and change with their 
influx and thus with the selected boundary condition. 
Furthermore, the open problem of its stability Iq intro- 
duces another uncertainty about their final sizes. 
IV. CONCLUSIONS 

In this work we develop a morphological phase dia- 
gram showing that during the collision of two Barchans 
four situations can be observed: coalescence, breeding, 
budding and solitary waves (Fig. ^ a, b, c, d and e). 
At large scales the collision process depicted here could 
lead to a selection of a characteristic size of dunes in 
a Barchan field. However, we only considered perfectly 
aligned dunes. 

Calculations of very large dune fields are still difficult 
because of high computational costs. One way out would 
be to consider a simplified model containing the main 
features of dune movement and interaction. For that 
purpose one could use the collision rules obtained in this 
work on a larger scale [23 • 



On the other hand. Figure [3 shows the inverse of the 
time during which the emerging dune does not split. The 
continuity of the curve at the 'budding-solitary wave' 
transition suggest a second order transition. Near the 
transition point the leaving dune splits in a infinite time. 
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